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Abstract: The spatial pattern of meteorological factors cannot be accurately simulated by using 
observations from meteorological stations (OMS) that are distributed sparsely in complex terrain. It is 
expected that the spatial-temporal characteristics of drought in regions with complex terrain can be better 
represented by meteorological data with the high spatial-temporal resolution and accuracy. In this study, 
Standard Precipitation Evapotranspiration Index (SPEI) calculated with meteorological factors extracted 
from ITPCAS (China Meteorological Forcing Dataset produced by the Institute of Tibetan Plateau 
Research, Chinese Academy of Sciences) was applied to identify the spatial-temporal characteristics of 
drought in Shaanxi Province of China, during the period of 1979-2016. Drought areas detected by SPEI 
calculated with data from ITPCAS (SPEL-ITPCAS) on the seasonal scale were validated by historical 
drought records from the Chinese Meteorological Disaster Canon-Shaanxi, and compared with drought 
areas detected by SPEI calculated with data from OMS (SPEI-OMS). Drought intensity, trend and temporal 
ranges for mutations of SPEI-ITPCAS were analyzed by using the cumulative drought intensity (CDI) index 
and the Mann-Kendall test. The results indicated that drought areas detected from SPEI-ITPCAS were 
closer to the historical drought records than those detected from SPEI-OMS. Severe and exceptional 
drought events with SPEI-ITPCAS lower than —1.0 occurred most frequently in summer, followed by 
spring. There was a general drying trend in spring and summer in Shaanxi Province and a significant wetting 
trend in autumn and winter in northern Shaanxi Province. On seasonal and annual scales, the regional and 
temporal ranges for mutations of SPEI-ITPCAS were different and most mutations occurred before the 
year 1990 in most regions of Shaanxi Province. The results reflect the response of different regions of 
Shaanxi Province to climate change, which will help to manage regional water resources. 
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1 Introduction 


Drought is a damaging and frequent natural disaster (He et al., 2011) that is significantly related to 
temperature abnormalities and shortage of precipitation (Salehnia et al., 20172). It is sensitive to 
climate change (Kogan et al., 2013; Jia et al., 2016) and has a significant impact on economic, 
ecological and agricultural activities. Research on the spatial-temporal characteristics of drought 
will help to identify the correlation of drought with meteorological, hydrological and ecological 
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processes (Yu et al., 2015; Yuan et al., 2015), and to further evaluate the impact of climate change 
on human society. 

Currently, research on the spatial-temporal characteristics of regional drought faces the following 
challenges: (1) the applicability of drought indices varies greatly due to the specific climatic 
features of the study area (Qiao et al., 2012; Yang et al., 2017); (2) the accuracy of drought 
monitoring depends heavily on the accuracy and spatial-temporal resolution of the input data for 
the drought index (Liu et al., 2016; Wang et al., 2016; Salehnia et al., 2017b); and (3) the occurrence 
and spatial-temporal distribution of drought events identified by the drought index need validation 
(Liu et al., 2018; Wang et al., 2018). It is imperative to find an appropriate drought index, a reliable 
validation method and the meteorological data with high spatial-temporal resolution to improve the 
quality of drought monitoring over specific regions. 

Meteorological drought is generally monitored and measured by multifactorial drought indices 
involving water and energy balance (Vasiliades and Loukas, 2009; Dai, 2011; Paulo et al., 2012; 
Gobena and Gan, 2013), such as the Palmer Drought Severity Index (PDSJ), Self-calibrated PDSI 
(Sc-PDSI) (Wells et al., 2004), and Standard Precipitation Evapotranspiration Index (SPEI) 
(Vicente-Serrano et al., 2010a). PDSI requires parameters including soil water holding capacity, 
groundwater level and lake/reservoir level (Ma et al., 2015), which are not available in many 
circumstances. Furthermore, the temporal scale of PDSI is fixed (Vicente-Serrano et al., 2011), 
although drought is a multi-scale phenomenon (Changnon and Easterling, 1989; Pandey and 
Rafmasastri, 2010). Sc-PDSI is an improved version of PDSI (Wells et al., 2004; Liu et al., 2015), 
but the temporal scale is still not alterable. In contrast, SPEI combines the sensitivity of PDSI to 
evaporation changes and the multi-temporal nature of the Standardized Precipitation Index (SPI) 
(Vicente-Serrano et al., 2010a). Besides, it needs fewer parameters compared with PDSI when the 
Thornthwaite (Th) method is used to calculate evaporation. Thus, SPEI has been widely used to 
elaborate the spatial-temporal characteristics of regional drought and to reconstruct the spatial- 
temporal distribution of drought (Allen et al., 2011; Meza, 2013). 

Much research on drought quantification, monitoring and analysis has been conducted by machine 
learning models or spatial interpolation methods with metrological factors extracted from 
observations from meteorological stations (OMS) (Jiang et al., 2015; Rhee and Im, 2017; Wang et 
al., 2019; Zhang et al., 2019). However, the spatial interpolation of metrological factors (e.g., 
precipitation) brings significant errors due to the sparse observations and the complex terrain. 
Accuracy of SPEI forecasting by machine learning models with multisource data needs further 
theoretical investigation, particularly because SPEI may vary under different climatic conditions 
(Rhee and Im, 2017). It is more reliable to adopt meteorological datasets with high spatial-temporal 
resolution to directly quantify the drought based on drought indices (Vicente-Serrano et al., 2010b; 
Hannaford et al., 2011; Wang et al., 2016). Many studies have tested the applicability of various 
meteorological datasets in drought analysis at multiple scales under different climatic conditions. 
For example, Salehnia et al. (2017b) calculated eight drought indices with AgMERRA (Agricultural 
Model Intercomparison and Improvement Project based on the NASA Modern-Era Retrospective 
Analysis for Research and Applications) precipitation data and OMS data for drought monitoring. 
Results showed that the correlation between drought events and the former was higher than that 
between drought events and the latter. Naumann et al. (2014) investigated five different datasets and 
three drought indicators with regard to their capability to improve drought monitoring in Africa. 
They found that the extents of drought areas identified by drought indicators were different. Spinoni 
et al. (2019) constructed a dataset of meteorological drought events from 1951 to 2016 based on 
Global Precipitation Climatology Centre (GPCC) and Climatic Research Unit (CRU) time-series 
datasets with a spatial resolution of 0.5?. It can be concluded that the spatial resolution and accuracy 
of the gridded meteorological data and a suitable drought index were essential to identify the 
occurrence, intensity and spatial extent of drought. 

Multiple metrics (e.g., duration, intensity and area) of drought events were extracted from 
historical drought records and served as "ground truth" to validate the estimates obtained from 
drought indices (Spinoni et al., 2019; Wang et al., 2019). Drought detected from data with the high 
spatial resolution and accuracy can be better validated by drought records. In addition, 
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meteorological data with a long temporal range can help to depict the trend of drought. The temporal 
range (1979-2016) and spatial resolution (approximately 10 kmx10 km) of ITPCAS (China 
Meteorological Forcing Dataset produced by the Institute of Tibetan Plateau Research, Chinese 
Academy of Sciences) (Chen et al., 2011) can better meet those requirements of drought analysis 
compared with other commonly used meteorological products such as TRMM (Tropical Rainfall 
Measuring Mission) (0.25?x0.25?) and CMORPH (Climate Prediction Center Morphing Technique) 
(starting from 1998). Numerous studies have demonstrated the efficiency of ITPCAS in hydrological 
and meteorological applications across China (Guo and Wang, 2013; Xue et al., 2013; Wang et al., 
2018). 

In this study, the spatial-temporal characteristics of drought in Shaanxi Province of China were 
reconstructed by SPEI and ITPCAS data, and validated by historical drought records. Firstly, the 
SPEI of all grids with the spatial resolution of 10 kmx10 km in Shaanxi Province was calculated 
with the temperature and precipitation data of ITPCAS (termed as SPEI-ITPCAS) from 1979 to 
2016. Then, drought areas detected from SPEI-ITPCAS were validated by drought records 
extracted from the Chinese Meteorological Disaster Canon-Shaanxi, and compared with those 
calculated by OMS and ordinary Kriging interpolation (OKI). Finally, drought intensity, trends and 
mutations of SPEI-ITPCAS in Shaanxi Province from 1979 to 2016 were analyzed using the 
cumulative drought intensity (CDI) index and the Mann-Kendall (MK) test. This study aims to give 
more details about the spatial-temporal characteristics of drought in Shaanxi Province, especially 
in those ungauged areas which were neglected or poorly represented in previous studies (Jiang et 
al., 2015; Wang et al., 2019; Zhang et al., 2019). 


2 Study area and methods 


2.1 Study area 


Shaanxi Province is located on the boundary between the western arid zone and the eastern 
monsoon zone of China (31?42'-39?35'N, 105?29'-111?15'E; Fig. 1). Meteorological disaster, 
especially drought, frequently occurs with damaging effects on the ecological system and regional 
economy. Shaanxi Province can be divided into three regions with regard to geography and climate: 
the northern Loess Plateau with a cold, arid climate; the middle Guanzhong Plain with a mild 
climate; and the southern Qinba Mountainous region with a warm, wet climate (Liu et al., 2013). 
Qinling Mountains run across the middle of Shaanxi Province with significant climate discrepancy 
between the northern and southern sides of the mountains. The Loess Plateau in the northern part 
of Shaanxi Province (briefly, northern Shaanxi) has a sub-arid monsoon climate with average 
annual precipitation varying greatly from 300 to 600 mm. It is prone to both droughts and 
waterlogging. The Guanzhong Plain in the middle of Shaanxi Province (briefly, Guanzhong) has a 
warm sub-humid monsoon climate with average annual precipitation between 550 and 700 mm. 
The Qinba Mountainous region in the southern part of Shaanxi Province (briefly, southern Shaanxi) 
has a subtropical humid monsoon climate with average annual precipitation between 800 and 1300 
mm. Flash floods and geological disasters occur frequently in this southern region. The annual 
mean temperature of the entire Shaanxi Province is about 13°C, with regional means of 10°C in the 
north, 13°C in the middle and 15°C in the south (Wan et al., 2013). For administrative purposes, 
Shaanxi Province is divided into ten areas of prefecture-level city (APC). Yulin and Yan'an belong 
to northern Shaanxi; Baoji, Xianyang, Tongchuan, Weinan and Xi'an belong to Guanzhong; and 
Shangluo, Ankang and Hanzhong belong to southern Shaanxi (Fig. 1). 


2.2 Data sources 


Temperature and precipitation data (1979-2016) with the horizontal spatial resolution of 0.1? 
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Fig.1 Topography, administrative areas and distribution of standard meteorological stations of Shaanxi Province. 


N, northern Shaanxi (the Loess Plateau in the northern part of Shaanxi Province); G, Guanzhong (Guanzhong Plain 
in the middle of Shaanxi Province); S, southern Shaanxi (the Qinba Mountainous region in the southern part of 
Shaanxi Province). The two blue lines are the northern and southern boundaries of Guanzhong. 


were extracted from ITPCAS. This dataset was merged with meteorological data from the following 
sources: China Meteorological Administration data, reanalysis dataset of Princeton, Global Land 
Data Assimilation System (GLDAS) data, and Global Energy and Water Exchanges-Surface 
Radiation Budget (GEWEX-SRB) radiation data. 

The Chinese Meteorological Disaster Canon-Shaanxi briefly summarizes the start-end months 
or seasons, duration, severity, areas and financial loss of drought events. However, drought areas 
were recorded vaguely. For example, the spatial extent was recorded as parts of the APC, or even 
parts of northern Shaanxi, Guanzhong and southern Shaanxi. Furthermore, the severity of drought 
is hard to quantify from these descriptions. Drought events extracted from records of the Chinese 
Meteorological Disaster Canon-Shaanxi (Zhai, 2005) during the period of 1979-2016 are shown in 
Table 1. It can be seen that Guanzhong and northern Shaanxi were more frequently affected by 
drought from 1979 to 2016. 

In this study, metrological observations including precipitation and temperature during the 
period of 1979-2016 were downloaded from the China Meteorological Information Center 
(http://data.cma.cn/). Because these data were collected under rigorous quality control, no 
modifications were made for this study. 


2.3 Methods 


Drought characteristics of Shaanxi Province were analyzed by the following steps in this study. 
Firstly, SPEI was calculated with OMS at each standard meteorological station of Shaanxi Province 
and then interpolated into 10 kmx10 km grids by OKI (termed as SPEI-OMS). Then, SPEI was 
calculated by temperature and precipitation data extracted from ITPCAS at 10 kmx10 km grids of 
Shaanxi Province (termed as SPEI-ITPCAS). Both SPEI-OMS and SPEI-ITPCAS were calculated 
on seasonal and annual scales, and reclassified into drought and non-drought areas. The drought 
areas and period of drought occurrence were identified and then validated by comparison with 
historical records of drought events listed in Table 1. Secondly, spatial distributions of drought 
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intensity, trend and temporal ranges for mutations of SPEI-ITPCAS in Shaanxi Province during the 
period of 1979-2016 were analyzed on seasonal (SPEI 3) and annual (SPEI 12) scales using the 
CDI index and the MK test. 


Table 1 Historical records of typical drought events occurred in Shaanxi Province from the Chinese 
Meteorological Disaster Canon-Shaanxi 


Year Season Drought areas 

1991 Summer South of northern Shaanxi, Guanzhong and southern Shaanxi (especially in Shangluo) 
1994 Winter North of northern Shaanxi and Shangluo 

1995 Summer South of northern Shaanxi and Guanzhong 

1995 Autumn All regions" 

1997 Spring All regions" 

1997 Summer All regions" 

1998 Autumn All regions" 

1998 Winter All regions" 

1999 Summer All regions" 

2000 Summer Northern Shaanxi, northern Guanzhong and east of southern Shaanxi 
2001 Spring All regions" 

2001 Summer Northern Shaanxi, Guanzhong and east of southern Shaanxi 

2002 Summer All regions of Shaanxi Province except for Tongchuan and south of Yan'an 
2008 Winter All regions” 

2010 Autumn All regions” 

2011 Summer Part of northern Shaanxi and Guanzhong 


Note: " means the drought events occurred in all ten prefecture-level cities. Northern Shaanxi refers to the Loess Plateau in the northern 
part of Shaanxi Province; Guanzhong refers to Guanzhong Plain in the middle of Shaanxi Province; southern Shaanxi refers to the Qinba 
Mountainous region in the southern part of Shaanxi Province. 


2.3.1 SPEI model 


SPEI (Vicente-Serrano et al., 2010a) reflects the characteristics of drought or waterlogging by 
calculating the difference between precipitation and potential evapotranspiration (PET). Currently, 
PET is generally calculated by Penman-Monteith (PM) equation (Allen et al., 1994) or Th 
(Thornthwaite, 1937) parameterization schemes. The PM scheme needs unavailable parameters 
such as reflectivity of ground objects in each grid of Shaanxi Province, whereas the Th scheme 
needs only latitude and monthly temperature data (Tsakiris et al., 2007) which are available for this 
study. Thus, the Th scheme was selected in this study. The calculation of SPEI was described as 
follows: 

Firstly, monthly potential evapotranspiration (PET;) was estimated by the Th scheme using 
Equation 1: 


10T We 
PET, 2 16K (=) : (1) 


where PET; (mm) is the potential evapotranspiration of month i; K is a correction coefficient 
computed as a function of the latitude and month; T; (°C) is the mean temperature of month i; I is 
an annual heat index; and m is a polynomial with parameter I. 
Secondly, monthly water deficit or surplus D; was calculated by monthly precipitation P; and 
PET; 
D, =P - PET. (2) 
The calculated D; values can be aggregated at different time scales. The water deficit or surplus 


Di, in a given month i and year j at the time scale of k (k<12) months was defined as Equations 3 
and 4: 
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Dj- È Dra 5D, (i<k), (3) 
1=13-k+i l=1 
Di Y D, (izk), (4) 
l=j-k+1 


where D; ; (mm) is the water deficit or surplus in the month ! of year j; and Di, jı (mm) is the water 
deficit or surplus in the month l of year j-1. 

According to the Kolmogorov—Smirnov test, the log-logistic distribution was the best to fit SPEI 
at time scales from 1 to 24 months compared with Pearson III and Gamma distribution (Vicente- 
Serrano et al., 2010a; Yu et al., 2014). Thus, the log-logistic distribution was used to calculate SPEI 
after D; series were normally standardized. The probability distribution function f(x) and 
cumulative distribution function F(x) of the log-logistic distribution for Dj; series were defined in 
Equations 5 and 6, respectively: 


EN m um 
ott HE] ° 
F(x) = E 5) (6) 
x-y 


where a, f and y are the scale, shape and origin (y<D;<°°) parameters, respectively. 
Finally, SPEI can be calculated by Equations 7 and 8: 
C, - CW € CW? 
1-4 d, - dW? - dw?" 


SPEI-W 


(7) 


W =,/-2In(P) . (8) 


W used in Equations 7 and 8 is the standardized value of F(x). Constants used in Equation 7 are 
listed as follows: Co=2.515517, C,-0.802853, C2=0.010328, di=1.432788, d2=0.189269 and 
d3=0.001308. In Equation 8, if P<0.5, then P can be calculated as P=1—F(x), where P is the 
probability of exceeding the determination of Dj; whereas if P>0.5, thus P can be calculated as 
P=F(x). 

In this study, SPEI was calculated on seasonal (SPEI 3) and annual (SPEI 12) scales. Categories 
of drought and waterlogging classified by SPEI are shown in Table 2. 


Table 2 Categories of drought and waterlogging classified by SPEI (Standard Precipitation Evapotranspiration 
Index) 


Drought Normal Waterlogging 
Exceptional Severe Moderate Normal Moderate Severe Exceptional 
SPEI <-2.0 —2.0- -1.0 —1.0- -0.5 —0.5-0.5 0.5—1.0 1.0-2.0 20.0 


2.3.2 Ordinary Kriging interpolation (OKT) 

OKI makes an unbiased estimation of variables at target locations by a weighted linear combination 
of neighboring observations. It has been widely used as a popular Kriging technique to interpolate 
meteorological factors. Previous studies showed that OKI outperformed in precipitation 
interpolation compared with methods including Inverse Distance Weighted, Thiessen Polygons and 
Simple Kriging (Zhang and Srinivasan, 2009; Szczesniak and Piniewski, 2015). Details of using 
OKI for precipitation interpolation can be found in the original study of Goovaerts (2000). 

2.3.83 Cumulative drought intensity (CDI) 


CDI on seasonal and annual scales can integrate the drought severity and frequency of various 
forms of drought, providing the ways to evaluate drought over the long-term period (Kim et al., 
2002; Loukas et al., 2008). It can well indicate the cumulative impact of drought on human society 
in a specific region. 
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We calculated CDI according to the following steps: firstly, we calculated the SPEI-ITPCAS of 
each grid on seasonal and annual scales; secondly, we accumulated the drought (SPEI-ITPCAS 
lower than —0.5) according to the temporal range (Kim et al., 2002; Loukas et al., 2008). The larger 
the absolute cumulative SPEI-ITPCAS is, the greater the CDI is. Drought areas with cumulative 
SPELITPCAS less than —14.0 were calculated by summing the areas of the grids with high CDI 
values. 

2.3.4 Mann-Kendall (MK) test 


MK test (Mann, 1945; Kendall, 1975) was used in this study for trend and trend mutation tests. MK 
test is widely applied in trend analysis of meteorological and environmental factors, such as 
precipitation, temperature and runoff (Ahmad et al., 2015). The principle of MK trend test was 
described below. 

MK trend test assumed that the original time series xi-(x1, X2, ..., Xn) was random and independent 
with no significant trend. The statistical variable S was defined in Equation 9: 


s- S ¥ smn(x,- x). (9) 


i=1 j=i+1 
where sgn is a symbolic function; and x; and x;are the values of variables of the time series. If n=10, 
S can be approximately normally distributed. In addition, variable Z was used to show the results 
of significance test, which was defined as Equation 10: 


(S -1)/Jvar(S) (S » 0) 


Z= 0 (S=0). (10) 


(S+1)/,/var(S) (S <0) 


In a two-sided hypothesis test, positive and negative values of Z corresponded to the upward and 
downward trends, respectively. If |Z|<Z1-«/2 (a is the significance level), the original hypothesis was 
accepted, namely, the time series data had no significant trend; otherwise, the hypothesis was 
rejected, and the time series data had a significant trend. Moreover, If |Z|=1.28, 1.64 and 2.32, the 
significance test was accepted at significance levels of 9096, 9596 and 9996, respectively. 

In this study, we conducted the mutation test of MK according to the following descriptions. 

The MK trend variation test assumed that the original time series x;=(X1, X2, ..., Xn) was random 
and independent with no significant trend. The statistical variable UF; was the forward sequence 
and was defined as Equation 11: 


UF, Soa) k-12,.,n (11) 
S. Sy Ko) TI (12) 
i=l 
T +1 I (13) 
i 0 [eom ) 
g(s)- 56-0. k=2, 3... n, (14) 
var(S,) = desde k=2, 3 o% Ü (15) 


where S; represents the total numbers of x; preceding all x; and is defined in Equations 12 and 13; 
E(Sx) represents the expectation of Sx; var(Sx) represents the variance of Sx; and k represents the 
sequence number. For each xi, rirepresents the total numbers of x; which is smaller than x; (j<i). 
Statistical parameter UF; obeys the normal distribution. If |UF;|>U./2, the time series data had a 
significant decreasing or increasing trend. UF"; was calculated by the same Equations 11-14 with 
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the reversed times series of x; When the two curves UF; and UF’; intersected between the critical 
boundary (y=+U,), the moment of intersection was the trend variation begins. 


3 Results 


3.1 Validation of SPEI-ITPCAS and SPEI-OMS 


Figures 2 and 3 show the drought areas of Shaanxi Province detected from SPEI-ITPCAS and 
SPEI-OMS, respectively, for years and seasons corresponding to typical drought events. It can be 
seen that the distribution of drought areas detected from SPEI-ITPCAS and SPEI-OMS were 
similar in summer of 1991, spring of 1997, autumn of 1998, winter of 1998 and spring of 2001 
(Figs. 2 and 3), which accounted for 31.396 of the total seasons. In contrast, there were significant 
discrepancies of drought areas between the remaining seasons, indicating that meteorological data 
extracted from OMS and ITPCAS were not spatially consistent with drought monitoring for most 
of the seasons. As shown in Table 1, drought occurred in all 10 APCs of Shaanxi Province in nine 
seasons: springs of 1997 and 2001; summers of 1997 and 1999; autumns of 1995, 1998 and 2010; 
and winters of 1998 and 2008. In the above nine seasons, the distribution of drought areas identified 
from SPEI-ITPCAS (Fig. 2) was consistent with the historical records from the Chinese 
Meteorological Disaster Canon-Shaanxi (Table 1). In contrast, in summer of 1999 and winter of 
2008, drought areas identified from SPEI-OMS (Fig. 3) were different with the historical drought 
records from the Chinese Meteorological Disaster Canon-Shaanxi (Table 1), with the absence of 
north of Guanzhong. 

Table 3 lists the discrepancies of drought areas detected from different sources, i.e., historical 
drought records (Table 1) and SPEI-ITPCAS (Fig. 2), while Table 4 shows the discrepancies of 
drought areas between historical drought records (Table 1) and SPEI-OMS (Fig. 3). Compared 
with the seasons and regions in which drought events occurred, the discrepancies of drought areas 
between the two SPEI results and historical records happened in summer and in north of 
Guanzhong and northern Shaanxi frequently. Discrepancies of drought areas listed in Table 3 were 
similar to those in the same season listed in Table 4. However, Table 4 presents more discrepancies 
than Table 3. Discrepancies of drought areas shown in Table 4 were more likely to occur in Yulin, 
Tongchuan and Weinan, which were frequently affected by drought. In conclusion, the detection 
of drought areas from SPEI-ITPCAS presented a higher accuracy than that from SPEI-OMS. 
Therefore, SPEI-ITPCAS was adopted to analyze the spatial-temporal characteristics of drought 
in Shaanxi Province. 


3.2 Distribution of CDI 


Figure 4 shows that the cumulative SPEI-ITPCAS (SPEI-ITPCAS lower than —0.5) in Shaanxi 
Province from 1979 to 2016 was within —10.0 to —20.0 on seasonal and annual scales. On seasonal 
and annual scales, the cumulative SPEI-ITPCAS was within —10.0 to —14.0 in most regions of 
Shaanxi Province. When the cumulative SPEI-ITPCAS was less than —14.0, the region was defined 
as a typical area with high CDI values. On the winter scale, the cumulative SPEI-ITPCAS was within 
—16.0 to —20.0 in most regions of Shaanxi Province and the CDI values were higher than those on 
other seasonal and annual scales. Figure 4a indicates that areas with higher CDI values in spring 
were mainly distributed on the eastern and western boundaries of Guanzhong and in southern 
Shaanxi (A1 and A2), especially in southeast of Guanzhong and northwest of southern Shaanxi (A2) 
with an area of more than 2000 km? and high CDI values on seasonal and annual scales. In summer, 
areas with high CDI values were mainly distributed in Yulin (B1) and Shangluo (B2; Fig. 4b). In 
autumn, areas with high CDI values were located in Yan'an (C1) and Hanzhong (C2; Fig. 4c). On an 
annual scale, the CDI values in western Yulin (E1) and eastern Shangluo and Weinan (E2; Fig. 4e) 
were high. However, Figure 4d indicates that the CDI values in winter in the whole region were 
abnormally high, which is not consistent with the historical drought records (Zhai, 2005). Thus, the 
CDI in winter was not analyzed in this section. 
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Fig. 2 Spatial distributions of drought and non-drought areas in Shaanxi Province detected from SPEI-ITPCAS 
(Standard Precipitation Evapotranspiration Index calculated with data from China Meteorological Forcing Dataset 
produced by the Institute of Tibetan Plateau Research, Chinese Academy of Sciences). Drought areas were detected 
as SPEI-ITPCAS lower than —0.5 (Table 2). The two blue lines are the northern and southern boundaries of 
Guanzhong. 
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Fig. 3 Spatial distributions of drought and non-drought areas in Shaanxi Province detected from SPEI-OMS 
(Standard Precipitation Evapotranspiration Index calculated with data from observations of meteorological 
stations). Drought areas were detected as SPEI-OMS lower than —0.5 (Table 2). The two blue lines are the northern 
and southern boundaries of Guanzhong. 
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Table 3  Discrepancies of drought areas detected from historical drought records of Chinese Meteorological 
Disaster Canon-Shaanxi and SPEI-ITPCAS 


Year Season Discrepancy of drought areas 
1991 Summer North of northern Shaanxi 
1994 Winter East of Ankang 

2000 Summer North of Guanzhong 
2002 Summer Tongchuan and south of Yan'an 


Note: SPEI-ITPCAS, Standard Precipitation Evapotranspiration Index calculated with data from China Meteorological Forcing Dataset 
produced by the Institute of Tibetan Plateau Research, Chinese Academy of Sciences. 


Table 4 Discrepancies of drought areas detected from historical drought records of Chinese Meteorological 
Disaster Canon-Shaanxi and SPEI-OMS 


Year Season Discrepancy of drought areas 

1991 Summer North of northern Shaanxi 

1994 Winter Part of Weinan 

1995 Summer South of northern Shaanxi and north of Guanzhong 
1999 Summer North of Guanzhong 

2000 Summer North of Guanzhong 

2001 Summer Yulin 

2002 Summer Tongchuan, south of Yan'an and Yulin 

2008 Winter Tongchuan, Xianyang and Weinan 

2011 Summer Guanzhong and parts of northern Shaanxi 


Note: SPEI-OMS, Standard Precipitation Evapotranspiration Index calculated with data from observations of meteorological stations. 
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Fig.4 Spatial distributions of cumulative drought intensity (CDI) in spring (a), summer (b), autumn (c) and winter 
(d), and on an annual scale (e) in Shaanxi Province from 1979 to 2016. A1 and A2 represent regions with high CDI 
values in spring; B1 and B2 represent regions with high CDI values in summer; C1 and C2 represent regions with 
high CDI values in autumn; E1 and E2 represent regions with high CDI values on an annual scale. The two blue 
lines are the northern and southern boundaries of Guanzhong. 
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Frequencies of severe and exceptional drought in Shaanxi Province from 1979 to 2016 are shown 
in Table 5. The rank of drought frequency on different scales was as follows: 
summer>spring>annual>autumn; however, the discrepancies between different scales were not 
significant. From Figure 4 and Table 5 we can conclude that, drought occurred most frequently and 
seriously in summer, followed by spring. Southeast of Guanzhong and northwest of southern 
Shaanxi were the most sensitive areas to the impact of drought with the highest CDI value and the 
largest drought area in spring and summer and on an annual scale. In addition, eastern Yulin, 
southern Yan'an and Hanzhong were also likely to be affected in summer and autumn. 


Table 5 Frequencies of severe and exceptional drought (SPEI-ITPCAS lower than —1.0) and regions with high 
cumulative drought intensity (CDI) values on seasonal and annual scales in Shaanxi Province 


Season/annual Severe and exceptional . We. 

alk aronet (Ames) Regions with high CDI values 

Sprin 12.331 Part of Hanzhong (A1), with an area of about 500 km?; part of Shangluo, 

pring i Weinan, Ankang and Xi'an (A2), with an area of about 2000 km? 
oa 12.937 Part of Yulin (B1), with an area of about 820 km’; part of Shangluo, 

j Weinan and Xi'an (B2), with an area of about 780 km? 
Part of Yan'an (C1), with an area of about 1000 km?; part of Hanzhong 
Autumn ey (C2), with an area of about 700 km? 
i + 2. 
Aal 12,250 Part of Yulin (E1), with an area of about 440 km’; part of Shangluo (E2), 


with an area of about 490 km? 
Note: A1, A2, B1, B2, C1, C2, E1 and E2 are depicted in Figure 4. 


3.3 MK test of precipitation, temperature and SPEI-ITPCAS 


As shown in Figures 5 and 6, temperature presented a significant increasing trend (Z>1.28) in most 
regions of Shaanxi Province on seasonal and annual scales. Only part of Ankang showed an 
insignificant decreasing trend (—1.28<Z<0.00) in summer, autumn and winter. However, the 
significance and variation trend of precipitation and SPEI-ITPCAS in MK test varied in most 
regions of Shaanxi Province on the seasonal scale. Tables 6 and 7 showed the significance and 
variation trends of precipitation and SPEI-ITPCAS on seasonal and annual scales in Shaanxi 
Province, respectively. In spring, the upward trend of precipitation occurred in most regions of 
Shaanxi Province, except for parts of Guanzhong. Moreover, the temperature showed a significant 
upward trend at 9996 confidence level. The SPEI-ITPCAS showed a drying trend in most regions 
of Shaanxi Province, except for part of southern Shaanxi. The variation trends of SPEI-ITPCAS 
and precipitation were opposite in most regions, which indicated that the abnormal increase of 
temperature dominated meteorological processes, leading to a drying trend in most regions. 

In summer, the downward trend of precipitation together with the significant upward trend of 
temperature led to a general drying trend in most regions of Shaanxi Province. This drying trend 
was especially significant in parts of Yan'an and Guanzhong, and west of southern Shaanxi. Almost 
all of the five APCs in Guanzhong showed a significant drying trend. The situation was consistent 
with historical drought records, which indicated that drought in spring and summer was quite 
serious in Guanzhong and northern Shaanxi (Zhai, 2005). 

In autumn and winter, the significant upward trend of precipitation mainly occurred in northern 
Shaanxi, while the downward trend mainly occurred in Guanzhong and southern Shaanxi. The 
temperature showed a significant upward trend in all regions. The variation trends of SPEI-ITPCAS 
reflected a composite interaction of precipitation and temperature with remarkable spatial 
heterogeneity in Shaanxi Province. In conclusion, we found a significant wetting trend in northern 
Shaanxi and a significant drying trend in parts of Xi'an, Weinan, Hanzhong and Ankang. 

The variation trends of precipitation and temperature in northern Shaanxi, Guanzhong and 
southern Shaanxi on an annual scale were similar with those on autumn and winter scales. The 
wetting trend in northern Shaanxi was less significant than the upward trend of precipitation, and 
the drying trends in Guanzhong and southern Shaanxi were more significant than the downward 
trend of precipitation on an annual scale. The phenomena showed that the significant upward trend 
of temperature in Shaanxi Province counteracted the upward trend of precipitation (leading to a 
less significant wetting trend in northern Shaanxi) and interacted with the downward trend of 
precipitation (leading to a more significant drying trend in Guanzhong and southern Shaanxi). 
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Fig. 5 Spatial distributions of variable Z of Mann-Kendall (MK) test for precipitation, temperature and SPEI- 
ITPCAS in spring (al—a3), summer (b1—b3), autumn (c1-c3) and winter (d1—d3) in Shaanxi Province. The two 


blue lines are the northern and southern boundaries of Guanzhong. 
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Fig. 6 Spatial distributions of variable Z of MK test for precipitation (a), temperature (b) and SPEI-ITPCAS (c) 
on an annual scale in Shaanxi Province. The two blue lines are the northern and southern boundaries of Guanzhong. 


Table 6 Significance and trend of precipitation in Mann-Kendall (MK) test on seasonal and annual scales in 
Shaanxi Province 


Season/annual scale Z value Trend Regions 
Spring Z=0.00 Upward All of Shaanxi Province except for Xi'an and Weinan 
221.28 Upward" Parts of Yulin and Ankang 
Z<0.00 Downward Xi'an and Weinan 
ZS -1.28 Downward" Parts of Xi'an and Weinan 
Summer 220.00 Upward Tongchuan and parts of Yulin and Yan'an 
221.28 Upward" 


Parts of northern Shaanxi, Guanzhong (except for Tongchuan) and 


a Board southern Shaanxi 


ZS -1.28 Downward" Parts of Baoji, Hanzhong and Weinan 
All of Shaanxi Province except for parts of Weinan, Xi'an, 
Aname ZONO se Shangluo and Ankang 
Z>1.28 Upward” Yulin and Yan'an 
Z<0.00 Downward Parts of Weinan, Xi'an, Shangluo and Ankang 
Z< -1.28 Downward" Parts of Xi'an and Ankang 


All of Shaanxi Province except for parts of Weinan, Xi'an, Ankang 


Winter Z=0.00 Upward and Hanzhong 
Z21.28 Upward" Yulin, Yan'an, and parts of Tongchuan, Xianyang and Baoji 
Z<0.00 Downward Parts of Weinan, Xi'an, Ankang and Hanzhong 
Z< -1.28 Downward" Parts of Xi'an and Ankang 
Annual Z=0.00 Upward Northern Shaanxi, parts of Guanzhong and southern Shaanxi 
221.28 Upward" Parts of Yulin, Yan'an and Tongchuan 
Z0.00 Downward Parts of Guanzhong and southern Shaanxi 
Z< -1.28 Downward" Parts of Baoji, Xi'an and Weinan 


Note: * means that the confidence level of trend is above 90%. 


3.4 Temporal ranges for mutations of SPEI-ITPCAS on seasonal and annual scales 


As shown in Figure 7, the temporal ranges for mutations of SPEI-ITPCAS can be notably discerned 
on seasonal and annual scales in Shaanxi Province. For example, mutations of SPEI-ITPCAS 
occurred during the period of 1990-2000 in west of Shaanxi Province and east of Guanzhong in 
spring (Fig. 7a). In summer, they occurred during the period of 1990-2000 in west and middle of 
Guanzhong (Fig. 7b). In autumn, they appeared after the year 2000 in north of northern Shaanxi 
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Table7 Significance and trend of SPEI-ITPCAS in MK test on seasonal and annual scales in Shaanxi Province 


Seasomammuel Z value Trend Regions 
scale 

Spring Z=0.00 Wetting Parts of southern Shaanxi and Baoji 
Z21.28 Wetting" Parts of Ankang 
Zx0.00 Drying All of Shaanxi Province except for parts of southern Shaanxi and Baoji 
Z< —1.28 Drying” Part of Yulin, Yan'an, and middle and east of Guanzhong 

Summer Z=0.00 Wetting Parts of Tongchuan and south of Yan'an 
Z21.28 Wetting" 
«0.00 [cie All of Shaanxi Province EUR m pus of Tongchuan and south of 
Z< -1.28 Drying" West of southern Shaanxi, and parts of Yan'an and Guanzhong 

Autumn 2270.00 Wetting Northern Shaanxi, west and north of Guanzhong, and east of Hanzhong 
221.28 Wetting" Yulin and north of Yan'an 
Z<0.00 Drying Middle and eastern parts of Guanzhong and southern Shaanxi 
Z< -1.28 Drying” Parts of Xi'an and Ankang 

Winter Z=0.00 Wetting Northern Shaanxi, and parts of Gaunzhong and southern Shaanxi 
Z=1.28 Wetting” Northern Shaanxi, and parts of Gaunzhong and southern Shaanxi 
Zx0.00 Drying Parts of Xi'an, Weinan, Hanzhong and Ankang 
Z< -1.28 Drying" Parts of Xi'an, Weinan, Hanzhong and Ankang 

Annual 220.00 Wetting Yulin, and parts of Yan'an, Tongchuan and Xianyang 
Z21.28 Wetting" 
Z«0.00 Ding All of Shaanxi E s Mobius: parts of Yan'an, 
ZS -1.28 Drying" Parts of Guanzhong, Hanzhong and Ankang 


Note: * means that the confidence level of trend is above 90%. 
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Fig.7 Temporal ranges for mutations of SPEI-ITPCAS in spring (a), summer (b), autumn (c) and winter (d), and 
on an annual scale (e) in Shaanxi Province. The two blue lines are the northern and southern boundaries of Guanzhong. 
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(Fig. 6c). In winter, they occurred after the year 2000 in east and middle of Guanzhong and in parts 
of southern Shaanxi (Fig. 7d). On an annual scale, mutations of SPEI-ITPCAS occurred during the 
period of 1990—2000 in west and middle of Guanzhong, and after the year 2000 in parts of northern 
Shaanxi (Fig. 7e). In conclusion, the regional and temporal ranges for mutations of SPEI-ITPCAS 
differed on seasonal and annual scales, and most mutations happened before the year 1990. As 
Shaanxi Province has three types of terrain with comparably diverse climates (i.e., cold and arid 
plateau in the north, mild plain in the middle, and wet and warm mountainous regions in the south), 
the regional responses to climate change varied. The discrepancies of temporal ranges for mutations 
of SPELITPCAS indicated that the responses of the three sub-regions (northern Shaanxi, 
Guanzhong and southern Shaanxi) to climate change may be quite different. 


4 Discussion 


4.4 Meteorological data with high spatial-temporal resolution 


Validation of SPEI-OMS and SPEI-ITPCAS has proved that meteorological datasets with high 
spatial-temporal resolution can better reflect the characteristics of drought in regions with complex 
climate than OMS. Thus, comprehensive evaluations of the existing meteorological datasets in 
drought monitoring are necessary. Vicente-Serrano et al. (2010b) has proposed a worldwide drought 
index dataset based on SPEI calculated with 0.5? net CDF reanalysis data. However, it is too coarse 
for drought monitoring at the provincial scale. ITPCAS, TRMM and CMORPH-PDF (CMORPH 
calibrated by probability density function matching and optimal interpolation method) are 
commonly used for provincial drought analysis in China (Tian and Li, 2015; Lu et al., 2018). 
ITPCAS has a higher spatial resolution than TRMM and a longer time range than CMORPH. 
Moreover, numerous studies have shown the superior performance of ITPCAS in predicting 
precipitation in hydrological simulations in China (Guo and Wang, 2013; He et al., 2015). Kan et 
al. (2013) pointed out that the accuracy of spatial distribution of precipitation data derived from 
ITPCAS for western China was better than those of TRMM Multi-satellite Precipitation Analysis 
(TMPA3B42V6) and CMORPH. In the future, we will downscale ITPCAS with auxiliary 
environmental factors to finer resolution so that the spatial-temporal characteristics of drought can 
be better identified even more precisely. 


4.2 Assessment of SPEI 


According to the algorithms of SPEI, the accuracy of SPEI depends on the following three factors: 
(1) accuracy of precipitation data; (2) accuracy of calculated PET; and (3) option of probability 
density function to fit water deficit or surplus. Precipitation from ITPCAS has been proved to be 
effective in modelling applications (Wang et al., 2018). Many scholars have verified that the log- 
logistic function, which is a probability density function and is used for calculating the SPEI, fits 
well the distribution of water surplus (Beguería et al., 2014; Yang et al., 2016). Therefore, the 
algorithms used for calculating the PET, including PM equation and Th method, need further 
evaluation. The error of PET calculated by Th method has been shown to be high in winter (Duan 
et al., 2013; Parsons et al., 2019). Figure 4d showed that in Shaanxi Province, the CDI value on the 
winter scale was abnormally higher than those on other seasonal and annual scales. This anomaly 
happened because low temperature in winter caused errors in the calculation of PET by Th method. 
Thus, estimates of PET in winter calculated by Th method need to be adjusted. 


4.3 Adjustments of Th method in winter 


Figures 2 and 4d and Table 3 showed that drought areas in winter detected by SPEI-ITPCAS were 
relatively consistent with the historical drought records, while the CDI in winter was abnormal, 
indicating that Th method tended to overestimate the intensity of drought in winter. Tsakiris et al. 
(2007) proposed that the ratio of precipitation to potential evapotranspiration (briefly, P/PET) can 
be used as an appropriate drought index and as a way to distinguish arid and humid regions (Paulo 
et al., 2012). Because arid and humid regions are relatively stable, PET calculated by Th method 
in winter can be adjusted by P/PET values. 
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5 Conclusions 


In this study, temperature and precipitation data were extracted from ITPCAS to calculate SPEI, 
and the calculation results were then validated by historical drought records. Spatial-temporal 
characteristics of drought on seasonal and annual scales in Shaanxi Province during the period of 
1979-2016 were analyzed by SPEI-ITPCAS. Results showed that SPEI-ITPCAS was superior to 
SPEI-OMS in detecting drought areas of Shaanxi Province on seasonal and annual scales. The rank 
of drought frequency on different scales was as follows: summer>spring>annual>autumn. The CDI 
values in Shaanxi Province during the period of 1979-2016 were different on seasonal and annual 
scales. In spring and summer, there was a general drying trend in most regions of Shaanxi Province 
while a significant drying trend in Weinan and Xi'an. In autumn and winter, there was a significant 
wetting trend in northern Shaanxi while a drying trend in other regions. Mutations of SPEI-ITPCAS 
occurred before the year 1990 in most regions of Shaanxi Province on different time scales. 
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